December 5, 2008 11:52 Contemporary Physics GraphTopo 



Contemporary Physics 

Vol. 00, No. 00, June 2008, 1-20 



RESEARCH ARTICLE 



OO 

o 
o 

O 

Q 

in 



Manifestations of topological effects in graphene 

Jiannis K. Pachos a * 

(Received 00 Month 200x; final version received 00 Month 200x) 

Graphene is a monoatomic layer of graphite with Carbon atoms arranged in a two dimensional honeycomb 
lattice configuration. It has been known for more than sixty years that the electronic structure of graphene can 
be modelled by two-dimensional massless relativistic fermions. This property gives rise to numerous applica- 
tions, both in applied sciences and in theoretical physics. Electronic circuits made out of graphene could take 
advantage of its high electron mobility that is witnessed even at room temperature. In the theoretical domain 
the Dirac-like behavior of graphene can simulate high energy effects, such as the relativistic Klein paradox. 
Even more surprisingly, topological effects can be encoded in graphene such as the generation of vortices, charge 
fractionalization and the emergence of anyons. The impact of the topological effects on graphene's electronic 
properties can be elegantly described by the Atiyah-Singer index theorem. Here we present a pedagogical en- 
counter of this theorem and review its various applications to graphene. A direct consequence of the index 
theorem is charge fractionalization that is usually known from the fractional quantum Hall effect. The charge 
fractionalization gives rise to the exciting possibility of realizing graphene based anyons that unlike bosons or 
fermions exhibit fractional statistics. Besides being of theoretical interest, anyons are a strong candidate for 
performing error free quantum information processing. 
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It is a rare example in science to have theoretical developments and technological applications 
growing simultaneously at the same rate as with research in graphene [H, Q] • This two dimensional 
monatomic layer of graphite was theoretically first investigated by Wallace 0]. Experimentally it 
was successfully isolated at the University of Manchester [4J . Graphene consists of carbon atoms 
positioned at the vertices of a two dimensional honeycomb lattice. Its electronic properties can 
be modelled by free electrons that are allowed to tunnel from one site of the honeycomb lattice to 
a neighboring one. The simple geometrical configuration of the Carbon atoms gives graphene, in 
the microscopic level, a surprisingly rich collective behavior. The latter is witnessed by quantum 
effects so resilient that they can influence the behavior of graphene even at room temperature. 
Due to its unique character graphene can successfully simulate phenomena that are expected 
to appear in relativistic quantum theories, such as the Klein paradox [Bj], and it revealed new 
effects, such is the anomalous quantum Hall effect @, 0|. 

The desire to employ graphene for studying topological effects comes from the simplicity 
with which these effects are encoded in graphene. The Dirac-like description of graphene gives a 
natural starting point. By varying the geometry of graphene one can create topologically distinct 
configurations such as a sphere or a torus on which the effective Dirac equation is defined. 
Moreover, by perturbing the graphene lattice and the tunnelling couplings of the electrons one 
can give rise to effective gauge and scalar fields. These fields are naturally coupled to the Dirac 
fermions that model the behavior of graphene. One can encode vortices on the effective scalar 
field with vorticity determined by the flux of the gauge field going through the vortex, much like 
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Figure 1. The graphene molecule with a C atom placed at each site of a honeycomb lattice. The electron orbitals of each C 
atom are depicted. Hybridization provides three equivalent sp2 orbitals that create the strong a bonds with the neighboring 
atoms. This gives rise to the covalent "backbone" of graphene with honeycomb lattice geometry. The remaining 2p z orbital 
creates a half filled it bond. The latter provides free electrons that tunnel from site to site giving rise to the fascinating 
electronic properties of graphene. 

in superconductors. How the Dirac fermions respond to these topological defects is the central 
subject of this article. 

1.1 The structure of graphene 

Graphene is a large molecule of carbon atoms, C, that are strongly bound together in a flat 
configuration. The C atoms are positioned on the sites of a honeycomb lattice, as shown in 
Fig. [H Three of the four outer electrons of each carbon atom are employed to strongly bond 
with its three neighboring atoms via a bonds. The 2p z orbital of the fourth electron creates a ir 
bond with a neighboring C atoms. The a bonds provide the covalent "backbone" structure of 
graphene with honeycomb lattice geometry. Their strength is responsible for the flexibility and 
robustness of the lattice. The ir bonds give rise to graphene's unique electronic structure. Each 
7r bond is half filled allowing the p orbital electrons to tunnel from one atom to the neighboring 
one. Graphene should be considered as a many body system: several electrons are allowed to 
tunnel from site to site, while at the same time they satisfy the Pauli exclusion principle. In the 
following we use this simple model to study the properties of graphene. 

1.2 Topological properties 

Topology is a branch of mathematics that is concerned with the global properties of geometrical 
objects rather than with their local details. These properties are supposed to be unaffected 
by continuous deformations of the objects. Restricted in two dimensions we can consider a 
surface that can be compact, like a sphere or a torus, or it can be open, like a flat sheet or a 
cylinder. These characteristics are very general and they are invariant even if small geometrical 
deformations are caused on the surface such as ripples. Compact surfaces can be topologically 
equivalent to a sphere, or to a torus, if they have a hole, like a ring doughnut. More complicated 
topologies are possible, too. To distinguish them we assign an integer number called the genus, 
g, that gives the number of "holes" of each surface. Examples of compact surfaces with genus 
g = 0, g = 1 and g = 2 are seen in Fig. 0(a), (b) and (c), respectively. It is clear that g is 
invariant even under deformations of the surface geometry as long as no holes are created or 
destroyed. In that spirit a cup is topologically equivalent to a ring! 

It is a natural question to ask how graphene, a naturally appearing material, can possess such 
extraordinary properties as the ones related to topology. The main reason is that its chemical 
bonding is so strong that it induces a stability in the geometrical properties of the underlined 
lattice. Indeed, several stable geometric variants of graphene appear in nature. In all of them 
each carbon atom has precisely three neighbors. The lattice of these molecules mainly consists 
of hexagons with a number of additional pentagons or heptagons that induce curvature. For 
example, one can encounter cylinders (Fig. 0(d)), so called nanotubes, produced from a folded 
graphene sheet, with the number of open faces, N, being equal to two. Spherical configurations 
are also possible (Fig. EJa)) the so called fullerenes with N = 0. The same geometry is also 
met in the football, where twelve pentagons are present in between the hexagonal plaquettes 
so that a spherical configuration can be created. It is also possible to have crossing nanotubes 
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Figure 2. Two dimensional surfaces with different topologies characterized by the genus, g, or the number of open faces, 
N. The surfaces (a), (b) and (c) are compact with N = and g = 0, 1 and 2, respectively. The surfaces (e) and (d) are 
open with g = and N = 2 and 4, respectively. 

corresponding to a surface with four open faces, as seen in Fig. [5Je). Apart from the robust 
geometrical structure of the underlined lattice, the presence of long range quantum coherences 
are also responsible for the emergence of the topological character in graphene. Indeed, at low 
enough temperatures one can assume that graphene's electrons can travel coherently through 
the whole molecule, thus detecting its geometry or topology. This makes graphene the natural 
laboratory to study the interplay between topology and quantum physics. 



1.3 The index theorem 

Classical as they may be, topological configurations can have a dramatic effect on the quantum 
behavior of a system. In particular, they determine the possible quantum states a system can 
have. This striking connection is demonstrated by the index theorem, initially introduced by 
Atiyah and Singer in the sixties [§J . The index theorem provides a direct way to determine certain 
properties of the spectrum of general Hamiltonians. Its application is rather straightforward 
providing a good alternative to the direct diagonalization of the Hamiltonian. This can be a 
tantalizing task for large systems. It is, therefore, not hard to imagine that the index theorem 
can have a dramatic impact on theoretical and applied sciences [9|, [l0|] . 




Figure 3. The energy levels of a fermionic system are considered here to be symmetric around E = 0. For half filling systems 
there are half as many fermions as there are possible states. The lowest energy is obtained by occupying the negative energy 
states with a single fermion due to the Pauli exclusion principle. The application of a properly tuned potential difference, 
V, can cause the system to change its state, if there are zero modes to be occupied. This will give rise to an electric current. 
In the absence of zero modes no change will be observed and the system will appear as an insulator. 



In particular, the index theorem provides information about the number of zero energy states 
that are present in systems such as graphene and its geometric variants Hi, LL2[]. They are 



important for determining conductivity properties. To understand why let us look closer at 
the nature of these molecules. It is common that the employed Hamiltonians have symmetric 
spectra around zero energy, called the particle-hole symmetry. The negative part of the spectra 
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is called valence band and the positive part conducting band. As mentioned above, graphene has 
approximately half as many free electrons as the positions they can occupy. Due to the Pauli 
exclusion principle these electrons fill up the valence band states one by one until the zero energy 
point. This point corresponds to half of the possible occupation positions. When adding a small 
potential difference one can witness conducting or insulating behavior depending on the presence 
of zero modes. If there are zero modes then they can be occupied in response to the electric field. 
This change in the electronic states appears macroscopicaily as a current. If there are no zero 
modes the small shift in the energy will cause no change in the behavior of the system, thus 
appearing as an insulator. 

Apart from the conducting properties, the presence of zero modes can give rise to exotic char- 
acteristics. It can be shown by general arguments that vortices encoded on graphene can cause 
the fractionalization of charge. This collective behavior is known from the fractional quantum 
Hall effect. Drawing the analogy further, the vortices in graphene could actually be possessing 
fractional anyonic statistics, which is neither similar to bosons nor fermions [l3j]. It is an excit- 
ing possibility to have a physical system that needs mainly isolation from the environment to 
reveal such fascinating properties. Experiments that reveal such aspects of topological behavior 
are usually hard to perform due to the complexity needed for topological properties to arise. 
Graphene, however, could achieve this with minimal requirements. 

2 Elementary properties of graphene 

As we have seen, graphene is a molecule consisting of C atoms arranged on a two dimensional 
honeycomb lattice. The elementary plaquette of the lattice is a hexagon and the atoms are 
placed on the sites of the lattice. Graphene's electronic properties can efficiently be described by 
the tight binding model, where spinless electrons tunnel from site to site along the links of the 
lattice without interacting with each other. This simple model provides a surprisingly accurate 
description of graphene and gives rise to a wealth of phenomena. The tight binding Hamiltonian 
is given by 



where the operator a\ creates an electron at site i, aj annihilates it at j and the coupling J gives 
the strength of the tunnelling transition. To determine the spectrum of Hamiltonian (JTJ) we take 
advantage of its periodicity. For that we notice that the system repeats itself with respect to 
a unit cell that comprises of two neighboring sites, e.g. along the V3 vector, as seen in Fig. [H 
This periodic structure allows us to perform a Fourier transformation that reduces the energy 
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Figure 4. Graphene can be modelled by a honeycomb lattice, where electrons are tunnelling from site to site. This lattice 
has hexagons as the elementary plaquettes. It can be split into two triangular sublattices, indicated by black and white 
circles and denoted here by A and B, respectively. The doted rectangular denotes the unit cell of the lattice. 
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eigenvalue problem to the diagonalization of the Hamiltonian in a single unit cell. Indeed, for 
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E-jjj A(.i(k).4(k)) ( E 3j_„. ") ($>) , (2) 

where a^k) and ae(k) correspond to the Fourier transformed operators that are positioned 
on the left and on the right of the unit cell, respectively. Now, it is straightforward to find the 
eigenvalue of the energy for electrons with certain momentum. This only requires diagonalizing 
the two by two matrix inside the integral and yields 



E(k) = ± jy 1 + 4 cos 2 ^ + 4 cos ^ cos ^1 . (3) 

Common as they might look, these energy eigenvalues possess a unique property. They become 
zero for isolated points of momentum. These are called Fermi points and are positioned at 
K ± = 2tt/3(1,1/\/3). The presence of Fermi points makes the low energy behavior of graphene 
so special. At half fermionic filling and at zero external potential graphene is at zero energy 
in the sense that the valence band is completely filled. In this case the behavior of graphene 
is dominated by the behavior of the Hamiltonian ([2]) near the Fermi points. To determine this 




Figure 5. The spectrum, E, of graphene near a single Fermi point as a function of the small momenta, p x and p y . The 
conical shape corresponds to a linear relation between the energy and the momentum that characterizes the Dirac equation. 
The positive spectrum corresponds to the conducting band and the negative to the valence band. 

behavior we substitute k = K± + p, expand the Hamiltonian in powers of |p| and take into 
account only the first order terms by assuming that |p| is small. This gives the equation 

( H + L ) *W = (V (<x °p)* ) *W = 7 • = E ^ & 

where p = — i/tV, er = (cr x , cr y , a z ) are the Pauli matrices and H± is the Hamiltonian kernel 
that corresponds to each Fermi point K±. Moreover, -y are the Dirac matrices given by 



and for simplicity we have set 3 J/2 = 1. Eqn. @ shows that, due to the presence of isolated Fermi 
points, the low energy behavior of graphene can be described by a relativistic Dirac operator, 
P = 7 " Pi with energy dispersion relations given in Fig. [5j The corresponding eigenvectors are 
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Figure 6. The introduction of curvature on a flat sheet of graphene by cutting a 7r/3 sector centered at the middle of a 
hexagon and gluing the opposite ends. This, finally, creates a cone that has positive curvature. The resulting lattice has a 
single pentagon at the apex of the cone, while the rest of the plaquettes remain hexagonal. 



four dimensional spinors, (| K+,4) , | K + B) , | K^A) , | K_5)) T , where the components are the 
wave functions of the unit cell elements A and B for each of the Fermi points K±. 

It is surprising that these effective Dirac fermions do not have any mass term. This fact is 
responsible for the high mobility properties of graphene. The presence of the Fermi points is 
responsible for the semi-metallic behavior of graphene. That is, it behaves as metal due to the 
vanishing gap between the valence and conducting bands, as seen in Fig. [5l but with charge 
carriers in small numbers due to the vanishing density of states at the Fermi level. Finally, there 
are emerging pseudo-spin degrees of freedom, where the Pauli matrices are acting, as we initially 
started with spinless electrons. 

Of course, we have not discovered a fundamental relativistic field. We have only shown that 
the collective behavior of graphene electrons can be described by the Dirac equation. Indeed, the 
obtained effective velocity of the electrons is three hundred times smaller than the speed of light, 
c. The relativistic description of graphene's electrons gives rise to numerous new phenomena that 
can be tested in the laboratory, like the Klein paradox [e| , the anomalous integer quantum Hall 
effect 0] and so on. 



2.1 Curvature effect 

The next theoretical surprise is related to the way curved graphene is described. As may be 
naturally expected the low energy description of a curved graphene surface is given by a Dirac 
equation that is defined on the corresponding curved surface. But this is not the full story. The 
curvature induces an effective gauge field that needs to be introduced in the Dirac equation. The 
latter is equivalent to an actual gauge field with magnetic field flux going through the molecule 
surface. This is a remarkably simple way to couple gauge fields to Dirac fermions. 

To understand the mechanism that gives rise to the gauge fields let us examine how one can 
insert curvature to graphene. The simplest way is by cutting a tt/3 sector from a graphene 
sheet and then gluing the opposite ends of the lattice, as seen in Fig. [6J This process adds a 
single pentagon at the apex of the created cone. All the other plaquettes remain hexagonal. This 
minimal geometrical distortion of the honeycomb lattice adds positive curvature. The curvature 
can be measured by circulating a tangential vector, T, around the apex, i.e. 

T-cfr = £. (6) 

But this is not yet the end of the story! The insertion of a single pentagon forces us to connect 
two sites that are of the same type, e.g. B for the case in Fig. [6J This is different compared 
to the flat graphene, where every site A has only B neighbors or via versa. This discontinuity 
has a dramatic effect on the corresponding spinor that describes the deformed lattice. Indeed, 
when the spinor, ^f, is transported around the apex by an angle 2ir it is forced at some point to 
jump from a site B to a site B instead of a site A. This reminds us of the effect that a magnetic 
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field has on the wave function of an electron moving on a closed path. A full circulation gives 
a phase factor to the electron wave function that is pro portional to the enclosed magnetic flux. 
This effect is known as the Aharonov-Bohm effect [141 ]. Quantum mechanically, the magnetic 
field introduces a discontinuity: the wave function of the electron should simultaneously describe 
the processes of the electron remaining static or spanning a closed trajectory. This is rectified 
by introducing a vector potential term in the Hamiltonian that makes the theory consistent. 
Following similar steps one needs to add a non-Abelian vector potential term, A, in the effective 
Hamiltonian description of the curved graphene that compensates the jump in the components 
of the spinor, ^. The circulation of A along a loop, C, around the apex is given by 

A-dr=V, (7) 

c 1 

where r v is the second Pauli matrix that mixes the K + and K_ components of the spinor. This 
fascinating emergence of the effective gauge field is simply caused by the geometric distortion 
of the graphene lattice. The manifestation of the gauge potential is witnessed by the properties 
of geometric variants of graphene, such as the spherical fullerenes, and it is in the core of the 
topological effects described here. 

The resulting Hamiltonian includes the curvature effect and a non-Abelian gauge field that 
couples the Fermi points, K + and K_. It is possible to rotate this Hamiltonian so that the gauge 
field reduces to two Abelian fields that act independently on each Fermi point. In this rotated 
frame the effective Hamiltonian is given by 

The covariant derivative is given by V = d — ifl, where fi is the spin connection that describes 
the effect of curvature on the vector character of the spinors. It has circulation 

n.£fr = -V, (9) 

c 6 

where a z acts on the A and B components of the spinor and C is a loop around the apex of a 
cone. The Abelian gauge potential satisfies 

A-dv = ^, (10) 
c 1 

with the corresponding magnetic field given by B = V x A. The zweibeins, e^, are introduced 
to make the Hamiltonian covariant in the induced curved space. If the curved space, with metric 
g^ v , has coordinates x^ and the local flat space, with metric rj kl , has coordinates £ fc then the 
zweibeins are defined by e£ = dx^/d^ and they satisfy g^ v = e^e"r] kl . Practically, the zweibeins 
transform the Pauli matrices from the flat to the curved space in order to be contracted with 
the corresponding covariant derivative. From the metric j™ one can define the corresponding 
Christoffel symbols, VL, = \g <Jp {d IJ ,gvp + d v gp P — d p gp U ), while the scalar curvature is given by 

R = sTR^p, where R^ pa = d a Y% - d p T^ + - (11) 

The final Dirac equation, = E^, where H is given in ([8]), describes the behavior of a Dirac 
spinor, coupled to a gauge field, A, on a surface with curvature R. This model approximates 
well the low energy behavior of the carbon molecules 15|, 1^, 17, 19], which corresponds to 



states with large wavelength. The lattice spacing appears to be small for these states and the 
continuous limit can be safely taken. To be precise, the continuous limit corresponds to the limit 
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(a) (b) (c) 




Figure 7. (a) The C60 fullerene molecule, where there are 12 pentagons necessary to produce the spherical configuration. 
Nanotube of armchair type (b) and of zigzag type (c). It is not necessary to introduce pentagons for producing nanotubcs 
as they have zero curvature. 



of infinitely large molecules, where states with arbitrarily large wavelength can be supported. 
Nevertheless, we will see that the continuous approximation describes graphene molecules well 
even for relatively small lattices. Hamiltonian ([8]) has all the necessary ingredients that are 
responsible for the topological character of graphene molecules. It is now of interest to see what 
specific form the Dirac equation takes for the various geometric deformations of graphene. 



3 Geometric variants of graphene 

Among the most common geometrical variants of graphene are the fullerene [13] and the nan- 



otubes 21[. Fullerenes, also known as buckyballs, are compact molecules with zero genus. As a 
representative, the C60 has sixty carbon atoms positioned in a spherical configuration, as seen in 
Fig. [7(a). In order to create the curvature of the sphere one has to introduce twelve pentagons on 
an initially flat sheet of hexagons. This is the case for all buckyballs of arbitrary size as we shall 
explicitly see in the next section. On the contrary nanotubes are just folded graphene sheets in 
a cylindrical configuration. Nanotubes have zero two dimensional curvature so there is no need 
for the presence of any pentagons. Due to their lattice structure there are many distinct ways 
in which a nanotube can be folded. Each one is characterized by different boundary conditions. 
Examples are the armchair (Fig. [7(b)) and the zigzag configuration (Fig. [7(c)). 

Of course, one can imagine more complex molecules with positive as well as negative cur- 
vature. Negative curvature is present, for example, at saddle points and can be introduced by 
heptagons. The latter introduce an effective magnetic field flux that is exactly opposite to the 
one corresponding to pentagons. By employing pentagons and heptagons, one can generate any 
desirable graphene surface. These include compact surfaces of arbitrary genus as well as open 
ones. In the following we consider only open surfaces that can be produced by a normal section 
of compact ones. An example is the nanotube that is produced from a torus by cutting it along 
its smallest circle. 

The Hamiltonian that describes the low energy spectrum of each of these molecules is given 
by a Dirac equation similar to the one in ([8]). Now the metric, g^ u , describes the curvature 
of the surface of each molecule. The gauge field, A, is responsible for the local magnetic flux 
through every polygonal deformation. Thus, by keeping track of the geometry of the surface of 
the molecule as well as its polygonal deformations one is able to construct the corresponding 
Dirac operator that describes its low energy limit. Tedious as it may seem there is a simple way to 
determine the total flux going through the molecule from its general topological characteristics. 
This is achieved with the aid of the Euler theorem, which is presented in the following section. 
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4 Euler theorem 



As we have seen, to establish the low energy description of graphene molecules we need to 
determine the total number of plaquette deformations. If we know the total number of pentagons 
and heptagons present in a molecule of certain genus, then we can determine the total flux that 
goes through the molecule. There is a simple way to relate the number of these deformations 
to the topology of the corresponding molecule. Indeed, the Euler theorem 22| gives a relation 
between the structural details of a polyhedral lattice and its total topological properties. There 
are numerous proofs of this theorem, the most pedagogical ones are based on a reduction of the 
polyhedral lattice to simpler ones without changing its topological properties. 




Figure 8. (a) A simple lattice corresponding to genus zero surface, where the vertices, V, edges, E, and faces, F, are 
depicted, (b) Two pentagonal and two heptagonal deformation can be introduced on graphene lattice without changing its 
fiat geometry. 



Consider a lattice that lies on a surface that we initially take to be compact of a certain genus, 
g. An example of such a lattice with genus g = is seen in Fig. 15(a). Let V, E and F be the 
number of vertices, edges and faces of the lattice, respectively. The Euler characteristic, x, is 
given by 

X = V-E + F = 2(l- 5 ), (12) 

where the second equation consists the Euler theorem. We now apply this theorem to the case 
of graphene molecules. Each graphene vertex has exactly three links. Let us assume that only 
pentagonal or heptagonal deformations are present. Denote by 715, n@ and 717 the total number 
of pentagons, hexagons and heptagons in the molecule. Then the total number of faces is given 
by the sum of the different polygons, F = + n$ + n-j. The total number of vertices is given 
by V = (5rt5 + 6ri6 + 7nj)/3 as each fc-gon has k vertices and each vertex is shared by three 
polygons. Similarly, the total number of edges is given by E = (5n 5 + 6hq + 7n 7 )/2 as each edge 
is shared by two polygons. Substituting this into Eqn. (|12j) we obtain 

n 5 - n 7 = 12(1 - g). (13) 

There are many characteristics that spring out from this result. If equal numbers of pentagons 
and heptagons are inserted then they do not change the topology of the surface as they cancel 
out. Indeed, on a flat graphene sheet one can add two pentagons and two heptagons without 
changing the curvature of the molecule away from these deformations, as seen in Fig. [5(b). This 
is consistent with the effective gauge field description, where pentagons and heptagons have 
opposite flux contributions. On the other hand, non-trivial topologies necessarily introduce an 
imbalance between pentagons and heptagons. The genus zero configurations need an excess of 
pentagons, while high genus surface need an excess of heptagons. Genus one surfaces do not 
need any pentagons or heptagons at all as they are equivalent to a flat sheet. 
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As particular examples we see that Eqn. (|13h reproduces the known case of a sphere with g = 
that gives x = 2. This corresponds to the fullerenes that have n.5 = 12 and n-j = 0. For the case 
of the torus we have g = 1 giving x = 0. Thus, no pentagons or heptagons are required. If one 
wants to consider the genus g = 2 surfaces then x = — 2 for which = and n 7 = 12. Fullerenes 
are compact surfaces, while the rest of the geometrical variants of graphene are usually defined 
on open surfaces. To apply our study to them we need to insert open boundary conditions to the 
Euler characteristic. This is simply achieved by considering the corresponding compact surface 
and introducing normal cuts to it. As an example, a torus with g = 1 can be cut along its small 
circle thus creating a cylinder with two open faces. The type of polygons between these two 
shapes is the same. So the generalization of the Euler theorem gives 

n 5 - n 7 = 12(1 -g)- 6N, (14) 

where N is the number of open faces. The surgical procedure presented above, employed to 
produce open faces gives two open faces when the genus is reduced by one. In the case of 
nanotubes the surface has g = and N = 2 that corresponds to n.5 = nj = 0. The case with 
g = and N = 4 corresponds to two crossing nanotubes. In general, arbitrary configurations 
of genus g and of open faces N can be considered. Then the imbalance between the polygonal 
deformations is determined by relation (|14p . It remains to see how the total flux associated with 
the total number of pentagons and heptagons can determine the low energy spectrum of the 
molecule. In the next section we see how the index theorem can provide such an elegant relation. 



5 The Atiyah-Singer index theorem 



The index theorem gives an insight into the structure of the spectrum of certain operators, 
like the Dirac operator, without having to diagonalize them [1, [jj, 10]. This information can be 
derived from general properties of the operators and the geometry of the space in which they 
are defined. In particular, the theorem is concerned with eigenstates that have zero eigenvalue, 
which in our case are the zero modes. Abstract as it might first sound it is not hard to follow 
a heuristic proof of the index theorem. Here we employ the approach based on the heat kernel 
expansion that provides the main elements, applicability and limitations of the theorem. The 
index theorem allows us to obtain a relation between the zero modes of the Dirac operator that 
corresponds to graphene and the total flux that goes through its surface. As the latter is related 
to the genus of the surface through the Euler characteristic we finally obtain, as promised, a 
relation between the zero modes and the topology of the surface. 
Our starting point is a Dirac operator of the form 




where P is an operator that maps a space V+ onto a space VL, while P f maps V- onto V+. If 
P is an n x m matrix, then pt is a m x n matrix and V+, VL is the space of n, m dimensional 
vectors, respectively. For the Hamiltonian (JSj) obtained in the previous section,^) is one of the two 
sub-blocks of the diagonal differing by the sign in front of the gauge field, A, P is a differential 
operator and V± is the space of complex functions. 

As we are interested in the zero modes o£p, i.e. the solutions of the equation Ipty = 0, we can 
define the number of distinctive eigenstates of P with zero eigenvalue as v+ and the ones of P^ 
as v_. To facilitate the bookkeeping we introduce the chirality operator 75 as 



(16) 
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Its eigenstates have eigenvalue ±1 if they belong to V±. To simplify our calculations we consider 
the operator fl 2 that has the same number of zero eigenstates as fl, but has the additional 
advantage of being diagonal 




It is possible to show that the operators P^P and PPt have the same non-zero eigenvalues. 
Indeed, if there is a state u such that PP^u = An, for eigenvalue A 7^ 0, then 

PP^u = \u P^P{P^u) = A(P f u), (18) 

which means that also the operator P^P has the same eigenvalue, A, corresponding to the 
eigenstate P^u. However, this is not necessarily the case for A = as P^u could be zero by itself. 
Let us now evaluate the following quantity 

Tr( 75 e-^ 2 ) = Tr(e-* ptp ) - Tr( e - tppt ) = £ e~ u + - ^e~ tA -, (19) 

A + A_ 

where A+ and A__ are the eigenvalues of the operators P^P and PP^ , respectively, and t is just an 
arbitrary parameter. The first step in (|19p is due to 75 acting on the exponential: it contributes 
a +1 to the eigenvectors of P'P, as they belong to V+, and a —1 to the eigenvectors of PP', as 
they belong to VI. In the last step we evaluated the trace as a sum over all the eigenvalues of 
the corresponding operators. As we have shown, every non-zero eigenvalue of Ptp corresponds 
to an eigenvalue of PP^. Thus, all terms with non-zero eigenvalues cancel out in pairs and only 
the zero eigenvalues of each operator remain, finally giving 

Tr( 75 e-^ 2 ) = ^ + -^. (20) 

The difference between the number of zero modes is in general undetermined. We call it the 
index of the fl operator, i.e. index($5) = v + — z>_. The above derivation provides us with an 
unexpected clue: due to the cancellation of the non-zero eigenvalue terms in (119p the expression 
Tr(75e _t ^ ) is actually t independent! 

What we just presented is merely a definition. Here, we are actually interested in finding 
out what is the value of index(fl). For that we shall emp loy an alternative way of calculating 
Tr(75e _ *^ ) ) called the heat kernel expansion method [23]. It states that for a two dimensional 
compact surface and for general / and D operators one has the expansion 

Ti(fe- t£) ) = ^- t Y / t k/2 a k (f,D), (21) 

fc>0 

where Tr denotes the trace of matrices as well as the integration with respect to spacial coordi- 
nates and ak are expansion coefficients. For / = 75 and D =p 2 we should obtain (j20j) . that is 
an expression which is t independent. To achieve this, the coefficients should be zero for all 
k except for k = 2 for which all the t contributions in (I21j) cancel out. The value of 0,2 can be 
found from the first order term in the t expansion of the exponential. For this calculation we 

first evaluate that Jj) 2 = -g^V^V u + {[j' M ,'y v ] F i^ ~ l R > where F ^ = d n A » ~ d » A n is the field 
strength and is a covariant derivative with respect to gauge and reparametrization transfor- 
mations. The magnetic field is given in terms of the field strength by = ^e kfJlU F llu . One can 
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now calculate that the non-zero expansion coefficient, is given by 

a 2 = Tr{ 75 (-^, 7 Hi^ - ±R)} = 2JJB ■ dS. (22) 

The integration in this equation runs over the whole surface. Combining the results from the 
two independent calculations of Tr^e - *^ ) we have 

indexQp) = v+ - v_ = -L J J B • dS. (23) 

This is the final formula that demonstrates the index theorem. It relates the total flux that 
goes out of the surface to the number of zero modes of the If) operator. The curvature, R, does 
not appear in the index expression as 75 is a traceless operator. The absence of curvature is 
a characteristic of two dimensional surfaces. Finally, it is worth noting that Eqn. (|23|) makes 
sense only if the right hand side is an integer. This is indeed the case when we consider compact 
surfaces. Then jj B • dS gives the total magnetic monopole charge inside the surface which 
takes discrete values due to Dirac's quantization condition of monopoles [24| . In the case of 
non-compact surfaces the theorem is applicable only when there is no flux going through the 
open faces. 



6 Spectra of graphene molecules 

In order to apply the index theorem to graphene and its geometrical variants we have to de- 
termine the contribution of the effective magnetic field to (|23p. As we have seen the effective 
magnetic field that appears in a certain graphene configuration, with minimally distorted lat- 
tices, is directly related to the genus of the graphene surface. In particular, we know that each 
pentagonal distortion in the honeycomb lattice gives rise to a specific circulation of the vec- 
tor potential, namely § c A • dr = ir/2, and similarly for heptagonal distortions. By employing 
Stokes's theorem we can relate the circulation of the gauge potential around a plaquette to the 
flux of the corresponding magnetic field. Indeed, we have 

A-dr= ff B-dS, (24) 

Cp J J Sp 

where C p is the looping trajectory for which the circulation of A is evaluated and S p is the 
corresponding area. Thus, the total flux that goes trough the surface of the molecule can be 
evaluated if we know the total number of pentagons and heptagons. As heptagons contribute 
the opposite flux than pentagons we are only interested in the difference of their numbers. 
This is exactly determined by the genus of the surface through the Euler characteristic, i.e. 
6x = n 5 — n 7 = 12(1 — g) — QN. The total flux of the gauge field that goes through the surface 
of the molecule is, hence, given by 



The total number of zero modes is the sum of the contributions from each sub-operator of H in 
(JH]). Consequently, we obtain the index of the Hamiltonian, H, that describes the molecule by 
adding these contributions. This finally gives 



index(#) = u + - i/_ = 6(1 - g) - 3N. 



(26) 
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This is our final result that relates the number of zero modes present in a certain graphene 
molecule to the topological characteristics of its surface. Eqn. (I26h gives us the minimum number 
of possible zero modes. The exact number is obtained from the index if either v + = or = 0. 
This actually holds for most of the cases when the index is non-zero. 

We can now compare our results with the known cases of graphene molecules. For example, a 
fullerene has g = and N = 0. Thus, it is expected to have six zero modes. This is indeed the 
case as the fullerenes have two triplet modes with energies very close to zero. For the case of a 
nanotube we have g = and N = 2. This gives rise to index zero, which is in agreement with 



261 ] . Nanotubes have either zero modes that 



previous theoretical and experimental results [2£ 
satisfy v+ = V- as happens for the armchair nanotubes (see Fig. E(b)) or no zero modes at all, 
as it is the case for the zig-zag nanotubes (see Fig. EJc)). The balance between the chiral zero 
modes is a consequence of the symmetry between the two directions along the nanotube that 
forces the two chiralities to be symmetric. Finally, the number of zero modes of more complicated 
molecules can be obtained from Eqn. (|26p . 




C60 C180 



C540 C1620 



Figure 9. The low lying energy levels for the fullerenes C60, C180, C540 and C1620 obtained by direct diagonalization of 
their lattice Hamiltonians. Two triplets are obtained that correspond to the six zero energy modes in the continuum limit. 
As the size of the fullerene increases the energies of the low lying modes tend rapidly close to zero. The rest of the energy 
levels are separated with an energy gap of order one above and below the depicted energies. The energies are in units of the 
tunnelling coupling, J. 



When comparing this to the spectra of real molecules one has to realize that small deviations 
are possible. Indeed, we have employed the continuous approximation in our derivation, while 
actual graphene molecules are defined on a lattice. Strictly speaking their lattice nature precludes 
them from having exact zero modes. Nevertheless, we expect to obtain low lying states that tend 
to zero modes when the molecule is taken to be larger and larger, thereby approximating the 
continuous limit better and better. This is clearly shown in Fig. [9l where the low lying spectra 
of the fullerenes C60, C180, C540 and C1620 are plotted. For each molecule two triplets are 
obtained corresponding to the six zero modes. This is as expected in the continuous limit with 
energies that approach rapidly zero as the size of fullerene increases. The rest of the spectrum 
is well separated from these low lying states. 

This is a surprising relation between the topology and the spectrum of graphene molecules. 
In particular, it suggests that the zero modes need certain conditions in order to exist. An 
important ingredient is the presence of magnetic flux that is effectively inserted in graphene 
molecules by geometrical deformations. By inducing curvature with minimal distortions of the 
underlying lattice in the form of pentagons and heptagons we were able to relate the number of 
these deformations to the general topological characteristics of the lattice surface. This resulted 
in Eqn. (|26p that relates the zero modes of a general graphene molecule with the genus and the 
number of open faces of its surface. 



7 The next step: Vortices 



The previously described topological effects are not the only ones that can be encountered in 
graphene. Recently, an inspiring relation between vortex-like topological defects and the number 
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of zero modes was demonstrated 0, 0, M, M, Ejll . [12] • This construction is not very different 
from the one we studied in the previous section. The crucial new element is the presence of a 
scalar field in the Dirac operator on which vortices can be encoded. This scalar field is generated 
in the effective picture by appropriate distortions of the tunnelling couplings of the graphene 
molecules. The main characteristic of the vortices is the presence of magnetic flux attached to 
them. The magnetic field associated with the vortices can be, in principle, introduced by external 
sources. Here we focus on the case where the field is generated by the geometrical deformations 
that we already introduced in the previous sections. This is a surprisingly simple mechanism 
compared to introducing external fields. 



Figure 10. The Kekule distortion introduced on the honeycomb lattice. A third of the hexagonal plaquettes has all the 
tunnelling couplings the same, while the rest two thirds of the hexagons have alternating strong and weak couplings. 



Let us first see how the scalar field can appear in the effective low energy behavior of graphene. 
Consider a particular distortion of the tunnelling couplings of graphene, namely the Kekule 
distortion, that is given by 

SJ T)k = i$( r ) e iK + .v, e i(K + -K_).r + c c ^ (2?) 
o 

where $ is a slowly varying function, K± are the Fermi points, r gives the position of the lattice 
sites and Vk, for k = 1,2, 3, indicate the positions of the three neighboring sites, as depicted in 
Fig. [H It creates a pattern of weak, J, and strong, J', couplings as represented in Fig. [TUJ The 
effect of the distortion in the Dirac equation is to introduce a complex field, in the following 

way 

where $ = \$\ e i2e f A ' dr . As we see, the scalar field is directly coupled to the gauge field. The 
vorticity of the scalar field is defined as 

if if - $d$* 

Vort($) = — * dr ■ dArg $ = — (h dr — -r= . (29) 



It was argued by Jackiw and Rossi [33f] and then proven by E. Weinberg [34|, [35|] that the index, 
which gives the zero modes of the Hamiltonian (|28p . is equal to the vorticity of the scalar field, 
i.e. index(.ff) = Vort( ( I ) ). An extension of this theorem to compact surfaces was given in [3(J. A 
direct calculation of the vorticity of the scalar field gives us 

mdex(H) = i jj B • dS. (30) 



Note, that this index is with respect to four dimensional Hamiltonians in contrast to the pair of 
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two Dirac equations we encountered in the previous. The index is now twice the previous one 
reflecting the fact that the scalar field has twice the charge of the fermions, 2e. The contour S 
is taken to enclose all vortices imprinted on the scalar field. 

Chemists have long surmised that in fullerenes not all the nearest-neighbor hopping couplings 
are equal. In numerical calculations it is found that the molecules can lower their electronic 
energy by undergoing a Kekule distortion. In particular, in the case of "leapfrog" fullerene 
molecules [36| C6o+6fc with k = 0, 1, ... it is possible to assign a Kekule distortion in a consistent 
way. Thus, in addition to the gauge field that corresponds to the pentagons one obtains a scalar 
field from the possible coupling distortion. 

From the previous study we know that the total effective magnetic flux of fullerenes is 3 • 2tt. 
This property remains unchanged when a Kekule distortion is introduced. The equality between 
vorticity (|29p and the corresponding index (|30p assign to the scalar field a total winding number 
of six. This is composed out of twelve half vortices attached to the twelve pentagons of the 
carbon lattice. The vorticity tt of each vortex, instead of the flux tt/2 of each pentagon, is due 
to the charge 2e that corresponds to the scalar field. 

The index theorem applied to fullerenes with coupling distortion seems different than the 
one we employed for the undistorted molecules. However, both give the same total number of 
zero modes. This underlines the fact that introducing a scalar field in a Dirac operator does not 
change its total index. The only extra complication that arises is the angular discontinuity of the 
scalar field around the vortices due to their ir vorticity. This forces us to introduce branch cuts, 
i.e. lines between pairs of vortices, where the scalar field takes zero value. Moreover, regularity 
of the vortices demands that the scalar field is smooth around the vortices. 




Figure 11. The coupling configuration of the C60 molecule, where vortices reside on the pentagons, (a) The Kekule dis- 
tortion, (b) Cuts between vortices are introduced by replacing a double bond with single ones, (c) An enlargement of the 
vortices is introduced by removing all double bonds connected to the pentagons. For C60 this removes all double bonds. 




Figure 12. The spectrum of the C60 molecule as a function of the double bond coupling, h, with single bond coupling equal 
to 1. There are six modes (two triplets) that are near zero energy, (a) With Kekule distortion, (b) With cuts between pairs 
of vortices, (c) With enlarged vortices, where all double bonds are removed. 



In Fig. [TlT a) we see the C60 fullerene with the Kekule distortion depicted in terms of single 
and double bonds. The vortices are placed at the center of the pentagon plaquettes. In Fig. [TiT b) 
branch cuts are introduced between connecting pairs of vortices. In Fig. [TTTc) we increased the 
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Figure 13. The spectrum of the C180 molecule, (a) With Kekule distortion, (b) With cuts between pairs of vortices, (c) 
With enlarged vortices. 



size of the vortices by removing the double bonds immediately connected to the pentagons. 
Fig. [12] shows the energy levels of the C60 fullerene as a function of the strength of the Kekule 
distortion for the coupling configurations depicted in Fig. 1111 In Fig. [13] the corresponding plots 
are presented for the C180 molecule. In both plots we see that two lines corresponding to two 
triplets have energies close to zero. These low lying states are interpreted as zero modes in the 
continuous limit. We clearly see that the insertion of branch cuts and the increasing of the size of 
the vortices makes the low lying energy states to be even closer to the zero energy. The difference 
of their energy from being exactly zero is due to the small size of the system. It is expected to 
converge to zero when larger molecules are employed 37]. 

To summarize, we showed that an effective scalar field is induced in the low energy description 
of graphene originating from a distortion in its tunnelling couplings. Importantly, the effective 
gauge field that originates from the insertion of pentagons and heptagons in the lattice couples 
directly to the scalar field imprinting vortices. As a result, the fullerene- type molecules have six 
half vortices. By employing the index theorem we demonstrated that there should be one zero 
mode for each pair of such vortices. While the effective gauge field is responsible for the presence 
of the zero modes the scalar field assures that the corresponding degenerate states are not locally 
distinguishable and remain localized around the vortex. 



8 Charge fractionalization 

The topological interplay between gauge and scalar fields in the relativistic description of 
graphene has, yet, another surprising consequence. It leads to the fractionalization of charge 
to the value e/2 for each 2ir vortex. Thereby, it gives a natural setting, where topological effects, 
such as anyons, can emerge. The charge fractionalization is not an actual mechanical breaking 
of the electrons. It is rather a collective effect, appearing in the quantum state of the system, 
which conspires to assign half a charge to each vortex. In the case of finite systems fractionalized 
charges, e/2, should come in pairs as the total charge should always be an integer. This is the 
case as, for a compact system, 2ir vortices came always in pairs. The fractionalization of charge is 
defined with respect to the difference between the electron density with and without the vortex. 
Nevertheless, it is not an uninteresting statistical effect that appears, e.g. due to averaging the 
wave functions that correspond to the presence and the absence of an electron. What actually 
happens is that the wave function of the system has fractional charge as an eigenvalue. 

In this Section we present how the coupling between electrons and a vortex, as described in 
the previous section for graphene, gives rise to charge 1/2. For that we employ the quantum 
field theory description [3a |. Consider the Dirac operator, H(Q V ), with a scalar field that has 
a full vortex. The complete set of orthonormal eigenstates of the Hamiltonian is then given by 
parameterized by the energy eigenvalue E. The index theorem states that the presence of 
vortices gives rise to zero modes. We assume that our system with a full vortex has only a single 
zero mode, *&o, related to it. We further assume that the eigenvalues E are symmetric around 
zero, i.e. the system possesses particle-hole symmetry. We also encountered this condition in 
Eqn. (|18h of the index theorem. 

In quantum field theory terms the state of the system is described by the quantum field 
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operator, ^. For the case of a vortex this operator can be expanded in the eigenmodes of the 
corresponding Dirac operator, ^ E , in the following way 

$ v = J2(pE*l + J E *l E ) + a*o. (31) 

E 

Due to the symmetry of the non-zero energy states, the modes ^> E and ^ V _ E can be paired. The 
first one appears with the particle annihilation operator, b E , that acts on the conducting band 
and the second appears with the anti-particle creation operator, d E , that acts on the valence 
band. However, the zero mode does not have a partner. It is represented simply by the fermionic 
operator a. The zero mode, as all the fermionic modes, can be either occupied or empty. Since 
both states have zero energy, the zero mode has a double degeneracy. Denoting the empty zero 
mode state by |— ) and the filled one by |+) we can represent the creation, a', and annihilation, 
a, operators as 

a\+) = |-), 0+1+) = 0, a\-) = 0, at]—) = |+). (32) 

In order to determine the total charge of the system we employ the charge operator, Q. This 
operator is given by 

Q = \ j d 2 x(i> ] i> -i>¥), (33) 

where the Schwinger's prescription is employed for its normalization. For the state (|31f) the 
charge operator becomes 

Q v = ^2(b E bE-d E d E )+a^a-^. (34) 

E 

Due to the balance between the energies in the conducting and valence bands, their contributions 
to the charge of half filled systems cancel out. Thus, the eigenvalues of the charge operator depend 
only on the state of the zero mode, giving 

Q v |-> = -^|->, Q v l+> = ^l+>- (35) 

This demonstrates that the vortex charge is 1/2 in units of the electron charge! Note that in the 
absence of the zero mode the neutrality condition, Q = would arise and no fractionalization 
would take place. 

The above result demonstrates that the fractional value of the charge is indeed an eigenvalue 
of the charge operator. Our derivation was simply based on the spectral symmetry of the Dirac 
operator and the presence of the zero mode in the vortex sector, which is topological in its origin. 
Hence, the result is completely general regardless of the shape of the vortex or other details of the 
system. Fractional quanta of charge have been experimentally measured in the case of fractional 
quantum Hall effect [H, H(J . 



9 Anyonic statistics 

The fractionalization of charge in two dimensions gives rise to another exotic effect. It makes it 
possible to realize anyons. The latter are quasiparticles, such as vortices, with fractional statistics 
that can range continuously from bosonic to fermionic. This is manifested by a non-trivial phase 
factor acquired by the state of two anyons when one is circulating the other, i.e. after two 
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successive exchanges, as seen in Fig. [T4j In the case of bosons or fermions this phase factor 
would be 1 either arising from bosonic, (+1) 2 = 1, or from fermionic, (— l) 2 = 1, statistics. 
In the case of anyons, the phase resulting from two successive exchanges can actually take any 




Figure 14. (a) Two quasiparticles are clockwise exchanged. The wave function of the quasiparticles is appropriately evolved 
according to their mutual statistics, (b) A succession of two exchanges gives rise to a single circulation of one particle around 
the other. In three spatial dimensions this circulation can be continuously transformed into a smaller loop and finally to a 
trivial one, where no circulation takes place. The latter evolution is equivalent to no transportation at all so the evolution 
operator should be the identity. Hence, the only possible exchange statistics is the bosonic, (+1) 2 = 1, or the fermionic, 
(— l) 2 = 1, as they are the only ones that give a trivial evolution after two exchanges. 



complex value. This possibility arises in two dimensions. There, one can think of anyons as 
composite objects of a magnetic flux and an electric charge. Then the phase factor acquired by 
the circulation can be attributed simply to the Aharonov-Bohm effect. Indeed, the charge of 
one of the anyons circulates the flux of the other one, as seen in Fig. [15j The acquired phase 
factor depends only on the number of circulations and not on other characteristics of the path 
as expected from a statistical evolution. 




Figure 15. In two dimension two exchanges lead to a trajectory that cannot be continuously deformed to the identity one 
due to the presence of the circulated particle. Thus, an arbitrary phase factor, e* s could be assigned to this evolution giving 
rise to anyonic statistics. Intuitively, this can arise from attaching an effective integer flux, ip, to each anyon and an effective 
fractional charge q. Then the Aharonov-Bohm effect is responsible for the phase factor e 1 " acquired by the wave function. 



Consider the case of graphene with Kekule distortion in the tunnelling couplings. This gives 
rise to a scalar field that can support vortices. Suppose that vortices can appear with vorticity 2ir. 
Then they will be accompanied by half a charge as we saw in the previous Section. Circulating 
such vortex around an identical one will give a phase factor n due to the Aharonov-Bohm effect, 



B • dS 



-2vr 
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This is a non-trivial phase that clearly demonstrates the anyonic nature of the vortices 4l| • 
The generation and detection of anyons in graphene is an exciting subject of current re- 

(I . If successful, it will give rise to a rich and versatile platform 



30, 31 



search [27], 128|, 

to perform experiments related to the fractional quantum Hall effect. Indeed, graphene offers 
several advantages, such as high electron mobility and resilience to temperature, that can prove 
valuable for probing these strongly correlated effects. Moreover, certain models with anyonic 
properties are of much interest in topological quantum computation. The latter promises to 
overcome the problem of environmental decoherence in the most efficient way. 



10 Conclusions 

In contrast to graphite, which is a rather soft substance, its constituent graphene layers have 
a robust lattice structure. This gives rise to graphene-based molecules that have a variety of 
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stable geometrical configurations, such as the fullerenes and the nanotubes. The predominant 
honeycomb lattice configuration of graphene is responsible for the effective relativistic behavior 
of its electrons. This behavior, combined with the long range quantum coherences present in 
graphene, allowed us to probe several exotic phenomena. Indeed, graphene and its geometrical 
variants provide an exciting laboratory for realizing numerous topological effects. 

Initially, we saw how the low energy of graphene can be described by the relativistic Dirac 
operator. An effective magnetic field naturally appears whenever we introduce a geometrical 
deformation on the honeycomb lattice. This magnetic field is naturally coupled to the Dirac 
fermions allowing the emergence of phenomena already known from the quantum Hall effect. 
By considering the underlined lattice we were able to relate local geometrical characteristics of 
graphene to the general topology of the molecule by employing the Euler character. As a result, 
the celebrated index theorem provides insight into the spectrum of graphene molecules just by 
considering their topology. This theorem is able to predict, in a simple way, the number of zero 
modes that are present in the geometrical variants of graphene. We compared these results with 
the spectra of fullerenes and nanotubes and found remarkable agreement. The spectra of more 
complicated molecules can be also predicted with our formalism. 

Subsequently, we considered the effect of tunnelling distortions on the relativistic properties 
of graphene. In particular, a scalar field emerges, on which vortices can be imprinted. An ef- 
fective magnetic field due to distortions of the lattice geometry can create non-trivial vorticity 
on the scalar field. We analyzed this effect for the case of fullerenes, where tunnelling distor- 
tions are expected to appear due to the simultaneous presence of pentagonal and hexagonal 
plaquettes. An appropriately modified index theorem can be employed that correctly predicts 
the number of zero modes. The number of zero modes is unchanged when the tunnelling dis- 
tortion is introduced. However, the wave functions of the zero modes are expected to become 
localized around vortices. This makes the exotic properties that accompany zero modes, such as 
the charge fractionalization, to become properties of the corresponding vortices. 

The generation of vortices can also appear in flat graphene sheets. There vortices with 2tt 
vorticity can be imprinted with the help of external magnetic fields. The resulting system behaves 
similarly to the fractional quantum Hall effect. The latter comprises of a continuum of two 
dimensional electrons at very low temperatures that are subject to a vertical magnetic field. 
There, vortices can support a variety of statistical behaviors. In the case of graphene we saw 
that the zero modes localized to vortices are responsible for the charge fractionalization as well 
as the emergent anyonic statistics of the vortices. It is an exciting perspective to realize such 
statistical phenomena with graphene. Its versatile behavior could prove valuable to technological 
applications of these many body effects. 
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